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ABSTRACT 

High-precision radial velocity surveys explore the population of low-mass exoplanets orbiting bright stars. This allows accurately 
deriving their orbital parameters such as their occurrence rate and the statistical distribution of their properties. Based on this, models 
of planetary formation and evolution can be constrained. The SOPHIE spectrograph has been continuously improved in past years, 
and thanks to an appropriate correction of systematic instrumental drift, it is now reaching 2 ms"i precision in radial velocity 
measurements on all timescales. As part of a dedicated radial velocity survey devoted to search for low-mass planets around a sample 
of 190 bright solar-type stars in the northern hemisphere, we report the detection of a warm Neptune with a minimum mass of 
16.1 + 2.7 M© orbiting the solar analog HD 164595 in 40 + 0.24 days . We also revised the parameters of the multiplanetary system 
around HD 190360. We discuss this new detection in the context of the upcoming space mission CHEOPS, which is devoted to a 
transit search of bright stars harboring known exoplanets. 
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1. Introduction 


Over the past 20 years, the radial velocity (RV) technique has 
benefited from continuous developments. Several major mile¬ 
stones in detection and characterization of extrasolar planets 
have been made possible with this technique, either in the 
context of large surveys or in follow-ups of transiting candi¬ 
dates. These developments allowed a few stabilized spectro¬ 
graphs like HARPS ( [Mayor et al.||2003| ) to gradually scan the 
area of Neptune-like objects and super-Earths and detect exo¬ 
planets with minimum masses as low as 1 Me ( jPumusque et al.j 


Recent exoplanet RV surveys show that over 40% of non¬ 
active solar-type stars probably harbor low-mass exoplanets 
( [Mayor et al.||2011t [Howard et al.[[2010l with a mass distribu¬ 
tion that seems to increase sharply toward low masses. In addi¬ 
tion, most of the lo w-mass exoplanets b elong to multiple com¬ 
pact systems (e.g., [Lovis et 'ar|[20111^ . These results overall 


^ Based on observations made with SOPHIE spectrograph on the 
1.93-m telescope at Observatoire de Haute-Provence (CNRS / OSU 
Pytheas), France (program 07A.PNP.CONS). 


agree with the observed properties of the transiting exoplanets 
detected by the Kepler Mission ( [Borucki & Koch|2011t[Howard[ 
[et al.[[201^ [Marcy et~ar[[2014[ ). Detection and characterization 
of new low-mass exoplanets orbiting bright stars is motivated by 
the need to obtain accurate measurements of the projected mass 
and the eccentricity and to identify all the components of multi¬ 
ple systems to constrain models of exoplanet evolution. Such ex¬ 
oplanets orbiting bright stars will also be key targets for a deeper 


and more extended characterization with CHEOPS ( [Broeg et^ 


[201 3[ [Eortier et al.[ 

2014J, for instance, which searches for the 

transit, and JWST ( 

Beichman et al.[[2014[), which searches for 


atmospheric signatures. 


The SOPHIE spectrograph ( [Perruchot et al.[[2008] ) on the 
1.93 m telescope of the Observatoire de Haute Provence was 
improved in June 2011 with the implementation of octagonal- 
section fibers ( [Perruchot et al.[[2QTT] ). This optimizati on led to 
an improvement in RV precision to nearly 2 ms“^ ( [Bouchy 


et al.[[2013] ) on timescales of a few tens of days. This preci¬ 


sion level is critical for detecting low-mass exoplanets around 
solar-type stars. In the context of the subprogram ''High preci¬ 
sion search for Neptunes and super-Earths'' of the SOPHIE ex- 
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trasolar planets research consortium ( [Bouchy et al.||2009] ), we 
redefined a sample of 190 bright stars on the following criteria 
: 1) in the northern hemisphere, 2) with spectral type G and K 
and 0.6 < B - V < lA, 3) on a volume limited to 35 pc, 4) 
non-active (vsini <4.5 kms“\ logR'^^ < -4.8), 5) not part of a 
binary system and not harboring giant planets at short orbital pe¬ 
riod, and 6) not part of RV su rveys carried out with HARPS-N 
GTO ( [Cosentino et al.|[MT^ at the TNG. This subprogram is 
intensively conducted since June of 2011. 

We here report on the detection of a Neptune-like exoplanet 
candidate orbiting the solar analog HD 164595. Section 2 de¬ 
scribes the SOPHIE observations and the data reduction. In Sect. 
3 we describe the correction of the systematic instrumental drift 
necessary to reach the 2 ms“^ precision level. This correction 
is then validated in Sect. 4 using the known planetary system 
HD 190360, for which we provide an updated orbital solution. 
Section 5 details our data analysis of HD 164595 and summa¬ 
rizes the orbital and physical properties of the planet. We finally 
present our conclusions in Sect. 6. 


2. Observations and data reduction 


Spectroscopic observations were conducte d with the SOPHIE 
spectrograph ( [Permchot et al.||2008 |2011| ). SOPHIE is a fiber- 
fed environmentally stabilized echelle spectrograph covering the 
visible range from 387 to 694 nm. We used the high-resolution 
(HR) mode {XIAX =75,000) with thorium-argon simultaneous 
wavelength calibration on fiber B. 

In total we gathered 75 measurements of HD 164595 with 
SOPHIE over 2.2 years. The exposure time was set to between 
800 and 1200 seconds to reach a signal-to-noise ratio (S/N) of 
140 at typical photon-noise RV uncertainties below 1.0 ms“^ 
while simultaneously averaging p-mode stellar oscillations. Tak¬ 
ing into account the wavelength calibration uncertainty (close to 
1.0 ms“^), the average RV uncertainty is 1.43 ms“^ Nine spec¬ 
tra were identified as low quality. One measurement has a S/N 
that is half of the average, indicating a strong atmospheric ab¬ 
sorption. Three measurements are affected by moon light con¬ 
tamination, which cannot be corrected with the simultaneous 
thorium mode. Five measurements present a flux anomaly of the 
thorium lamp, which prevented us from computing the simulta¬ 
neous drift. They were discarded for the data analysis in Sect. 
4.4. 

The s pectra were reduced and extracted using the SOPHIE 
pipeline ( [Bouchy et al.|[2009] ), and the resulting wavelength- 
calibrated 2D spectra were correlated using a numerical binary 
mask corresponding to spectral type G2 to obtain the radi al ve¬ 
locity measurement ( jBaranne et al.|1996{|Pepe et al.|2002| ). The 
full width at half maximum (FWHM), contrast, and bisector span 
of the cross-correlation function (CCF) are also provided by the 
pipeline. 

The complete radial velocity dataset, after the zero-point 
drift correction described in the next section, is reported in the 
Appendix and plotted in Fig.[^ 


3. Correction of the systematic instrumentai drift 

Despite a clear improvement in RV precision provided by the 
implementation of octogonal fibers on timescales of a few tens 
of days ( [Bouchy et al.||2013l [Permchot et al.[[20TT] ), we iden¬ 
tified long-term variations of the zero point of the instrument 
that added up to about + 10 ms“^ over the last 3.5 years (cf. 
Fig.[^ upper panel). The first ramp seems to be correlated with 
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the thorium-argon lamp aging. One jump of about -7 m/s is re¬ 
lated to the implementation of additional octagonal fibers after 
the double scrambling in December 2012 (BJD=55950). A sec¬ 
ond jump of -7 m/s at BJD=56775 appears after the coating of 
the secondary mirror of the 1.93 m telescope, which introduced 
a significant change in the flux balance across the spectral range. 
Other events, with a timescale of a few weeks, are clearly corre¬ 
lated with strong changes of the outside temperature, which are 
propagated to the tripod of the spectrograph at a level of a few 
tenths of degrees through the telescope pillar. A thermal control 
loop of the spectrograph tripod is under development to avoid 
such thermal conduction coming from the telescope pillar. 


To track this long-term zero-point drift of the instrument, we 
systematically monitored a set of RV constant stars each night, 
namely HD1851 44 ([Howard et al.|20TQl ), HD9407, HD221354, 
and HD89269A ( [Howard et al.[[2011[ ), which were monitored 
with HIRES on the 10.2m Keck telescope with an RV dispersion 
of 2.0, 1.7, 1.9, and 2.0 ms“^ , respectively, over several years. 
In addition to this set of four RV constant stars, we also used 
51 targets from our sample that have at least ten measurements. 
We recursively built an RV constant master using these two sets 
of stars. A first master was created using a running median on 
the four reliable constant stars. This master was then subtracted 
from the second group, and all the stars with a corrected RMS 
lower than a threshold of 3 m s“^ were included in the first sub¬ 
set. The offset for each star was adjusted as well to minimize 
the dispersion of the residuals. A double weight was given to the 
four reliable constant stars. This process was repeated until con¬ 
vergence, that is, when no additional star was added to the first 
subset. 


The running median was implemented with the criterium to 
average 15 measurements. This was preferred instead of taking 
a fixed time window to give the same weight to all correction 
points. Consequently, a period with a higher density of obser¬ 
vations will increase the temporal resolution of the correction. 
On the other hand, a period with sparse measurements will lead 
to a lower temporal resolution instead of a poorer constraint on 
the correction point. In practice, with this criteria of 15 measure¬ 
ments averaged, the typical timescale of the correction is nine 
days. 

The most satisfying correction uses a threshold of 3 ms“^ 
and converges after four iterations with a total of 23 constant 
stars in the final set. The radial velocities of these 23 constant 
stars as well as the derived RV constant master are displayed 
in the top panel of Fig. RV variations are clearly visible on 
different timescales, from a few days to a global trend over the 
last 3.5 years. This approach allowed us to gather enough data to 
correctly cover the time series and to average all physical effects 
such as stellar activity or signals of small planets, leaving only 
the instrumental systematics. 

To illustrate the gain provided by our correction of the 
systematic instrumental drift, the RV of the constant star 
HD 185144 is displayed in the bottom panel of Fig. For this 
target, the RV dispersion decreases from 5.4 to 1.6 ms“^ when 
the RV constant master is subtracted. Figure shows the RMS 
distribution of the 55 stars (including the four reliable constant 
stars) before and after correction of the RV constant master and 
illustrates the fact that the SOPHIE RV precision is close to 
2 ms“^ 


^ The RV constant master is available upon request. 
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Fig. 1. Upper panel : Zero-point drift correction (red line) as a function 
of time and the adjusted RV data set used to create it (blue dots). Bottom 
panel : RV time series of HD 185144 after the zero-point drift correction. 
This correction decreases the RMS from 5.37 ms“^ to 1.61 ms“^ . 



4. Effect on the planetary system HD190360 

4.1. Validation of the correction 

To validate the correction of the RV constant master, we ana¬ 
lyzed the known planetary system HD 190360 that was observed 
with SOPHIE in the context of the science validation of the 
octogonal fibers ( [Bouchy et al.|[MT^ . The plane tary system 
HD 190360 is composed of a long-period Jupiter ( jNaef et aLj 
|2003| ) and a warm Neptune ( jVogt et JT]|2005| ). A total of 41 
measurements of HD 190360 were gathered with SOPHIE us¬ 
ing the same setting as for HD 164595. The radial velocity table 
is reported in the Appendix. The RV dispersion is 9.5 m s“^ and 
9.3 ms“^ before and after correction of the RV constant mas¬ 
ter, respectively. The span of our observations was 2.4 years, 
which is about one-eighth of the period of HD 190360b. This 
does not allow efficiently constraining the orbital parameters of 
this planet. We therefore focus our analysis on the warm Neptune 
HD190360C. 

HD 190360 is a quiet star with a logR'^^ derived from SO¬ 
PHIE spectra of -5.13. No periodic signals in either the bisector, 
the CCE-EWHM, or the activity index are present at the two pe¬ 
riods found in RVs, which excludes stellar activity as the origin 
of the Keplerian signals. 


Fig. 2. Initial (blue) and final (red) distribution of the RMS in the 55-star 
sample. 


To study the impact of our zero point drift correction, we an¬ 
alyzed both the HD 190360 corrected and uncorrected data set 
with the same methodology. We used the planet analysis and 
small transit investigation software (PASTIS) ( jPiaz et al.|2014| ), 
which is a Bayesian analysis software, for fitting a two-planet 
model. HD 190360b orbit was not fixed but subjected to strong 
priors for the Bayesian analysis. Normal distributions were used, 
centered on the pu blished parameters w ith widths equal to the 
published errors in [Wright et al.j ( [20091 ). the other hand, the 
other planet parameters were let free, with a Jeffrey distribution 
for the period, a beta distribution for the eccentricity as defined in 
[Kipping[ ( [2013| ), and uniform distributions for the other parame¬ 
ters. We ran 20 MCMC chains with 300 000 steps each. Only the 
stationary parts of each chain were kept. The distributions of pa¬ 
rameter values of all the remaining uncorrelated chain links then 
correspond to the target joint posterior-probability distributions 
for the model parameters. 

Table reports our orbital parameters and their uncertain¬ 
ties for both cases. The results are v ery similar for the tw o data 
sets and agree well (within Icr) with [Wright et al.[ ( 2009[ ). How¬ 
ever, the errors in the uncorrected data are systematically higher 
(by up to 50%). The dispersion of the residuals is better in the 
corrected set with 2.1 ms“\ which is consistent with the pre¬ 
cision obtained on the constant stars, instead of 3.5 ms“^ with 
the uncorrected data. Moreover, the additional jitter necessary 
to set the reducedto 1 is 1.9 ms“\ instead of 2.9 ms“^ for 
the uncorrected data set. This test demonstrates the relevance of 
the zero-point drift correction, which significantly improves the 
quality of the RV data. 


4.2. Analysis of the joint HIRES and SOPHIE data 

Finally, we combined our data set with the Keck observations 
used by [Wright et al.[ ( [2009| ) to update orbital parameters for the 
HD190360 system. 100 MCMC chains of 300 000 steps were 
run with PASTIS to conduct a global search without a priori for 
a two-planet fit. The RV time series with our orbital solution is 
plotted in Fig. and the phase-folded radial velocities for both 
HD 190360 b and c are plotted in Fig.|^ The orbital and physical 
parameters of the two planets are reported in Table 

The updated orbital parameters all agree well (within 1 cr) 
or marginally well (within 2 cr) with the values previously pub- 
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Table 1. Comparison of the orbital solutions for HD 190360c with cor¬ 
rected and uncorrected SOPHIE observations. 


Parameter 

HD 190360c 

HD 190360c 


(uncorrected data) 

(corrected data) 

P [days] 

17.117 + 0.015 

17.127 + 0.011 

To [BJD] 

56149 + 4 

56134 + 4 

e 

0.088 + 0.1 

0 . 062 + 0 - 0^6 

0 ) [deg] 

237.1 +94 

248.5-*’ 

K[ms-'] 

5.72 + 0.79 

5.64 ± 0.55 

^meas 

41 

41 

o-(O-C) [ms-i] 

2.9 

2.1 

Jitter (0 - C) [ms“^] 

3.5 

1.9 



51000 52000 53000 54000 55000 56000 57000 


„ 0.010 - 
0.005 
0.000 
V -0.005 


. hv 


51000 52000 53000 54000 55000 56000 57000 

BJD - 2400000 [days] 


Fig. 3. Radial velocities of HD 190360 for HIRES (blue dots) and SO¬ 
PHIE (red dots) with the best-fit orbital solution (solid line). Error bars 
include the photon noise and jitter. 




lished in [Wright et al.| ( |200^ . The precision on all parameters is 
significantly improved except for the semi-amplitude and oj. As¬ 
suming a stellar mass of 0.96 + 0.1 M© as in [Vogt et al.| ( |2005) , 
the minimum mass of HD 190360 b and c is 1.50 ± 0.15 Mj^p and 
20.28 ±3.16 Me, respectively. 

The additional jitter (intrumental and stellar noise) neces¬ 
sary to set the reduced to 1 is 2.9 ms“^ for HIRES and 1.8 
m s“^ for SOPHIE. Additionally, the dispersion of the residu¬ 
als is significantly higher for the HIRES data (3.3 ms“^) than 
for SOPHIE (2.4 ms“^). To take into account the much longer 
time span of the Keck observations (11.9 years), we computed 
the median dispersion of the residuals for all the 2.4 year win¬ 
dows in the HIRES data. Its value is 3.1 m s“^ still significantly 
higher than for SOPHIE. This could be the result of either a 
higher instrumental noise or/and a higher stellar jitter at the time 
of the Keck observatio ns. The latter cas e seems unlikely since 
the logR^j^ reported in [Vogt et al.|(|2005|) is close to our value at 
-5.09. 


Fig. 4. Phase-folded radial velocities of HD 190360c {upper panel) and 
HD190360b {bottom panel) for HIRES (blue dots) and SOPHIE (red 
dots) with the best-fit orbital solution (solid line). Error bars include the 
photon noise and jitter. 


Table 2. Updated orbital and physical parameters of HD 190360 b and c 


Parameter 

HD 190360b 

HD190360C 

P [days] 

2867.9 + 7.7 

17.1186 ±0.0016 

To [BJD] 

59271 + 19 

55570 . 3+^5 

e 

0.343 + 0.017 

0.107 ±0.07 

OJ [deg] 

14.7 + 32 

305.8!|o 

K[ms-i] 

23.39 + 0.46 

5.20 ± 0.37 

msini [M®] 

475.16 + 49.0 

20.28 ±3.16 

N,„eas (SOPHIE) 


41 

N^eas (HIRES) 


107 

Jitter (SOPHIE) [ms”! 

] 

1.8 

Jitter (HIRES) [ms”'] 


2.9 

o-(O-C) (SOPHIE) [ms 

-'] 

2.4 

cr(O-C) (HIRES) [ms" 

■^] 

3.3 


5. HD164595 data analysis 

5.1. Stellar properties 

HD164595 is a G2V star with a magnitude V=7.08. Its Hip- 
parcos parallax (tt = 34.57 + 0.5 mas) implies a distance of 
28.93 + 0.4 pc. This star is considered to be one of the clos¬ 
est known solar analogs. Its parameters were derived by jPortoj 


de Mello et al.j ( |2014| ) in the context of their photometric and 


spectroscopic survey of solar twin stars within 50 parsecs of the 
Su n. From the SOPHI E CCFs and using the calibration given 
by jBoisse et al. ( 2010| ), we derived a vsini of 2.1 + 1 kms“^ 
HD 164595 is known to be a quiet star. [Wright et al.| ( |2004| ) re¬ 
ported a logR^j^ of -5.0 from observations made between 1998 
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Table 3. Observed and inferred stellar parameters for HD 164595. 


Parameter 
Spectral Type 
TT [mas] 
Distance [pc] 
V 
B-V 
[Mo] 
TeffUQ 
^Ogg 
[Fe/H] 

V sin i [kms“^] 
< > 


HD 164595 
G2V 

34.57 + 0.5 
28.93 + 0.4 
7.1 
0.635 

0.99 + 0.03 
5790 + 40 
4.44 + 0.05 
-0.04 + 0.08 
2.1 + 1 
-4.86 + 0.05 


Reference 

Hipparcos 

Hipparcos 

Hipparcos 

Hipparcos 


Porto de Mello et al. 

(2014 

Porto de Mello et al. 

(2014 

Porto de Mello et al. 

(2014 

Porto de Mello et al. 

(2014 

Porto de Mello et al. 

(2014 


this paper 
this paper 


and 2003. More recently, [Isaacson & Fischer] ( |2010| ) reported a 
slight trend of logR'^^ from -4.97 to -4.86 between 2005 and 
2009. Following the method established by Boisse et al. (2010), 
we derived from our SOPHIE spectra an average value of -4.86 
+ 0.05 for logR^j^ . Furthermore, computing the median value 
for each season, we find that the log slightly increased from 
-4.91 in 2012 to -4.88 in 2013 and finally -4.81 in 2014. We also 
note that the bisector span shows an increasing dispersion over 
three seasons from 1.7 ms“^ in 2012 to 2.0 ms“^ in 2013 and 
finally 3.3 ms“^ in 2014. During the last season, the bisector 
span is marginally correlated with the RV. All these factors seem 
to indicate that the stellar activity of HD 164595 evolves along 
a magnetic cycle toward increased activity. Table [^summarizes 
the stellar parameters of HD 164595. 


5.2. Distinguishing systematic signals and aliases 

The Generalized Lomb-Scargle periodogram of the uncorrected 
radial velocities of HD 164595 is displayed in the upper panel 
of Fig. A number of significant features appear, the most 
interesting being a triplet of peaks at 40+5 days. The high¬ 
est one lies below the 10“^ false-alarm probablity (FAP) level. 
The FAP was calculated with Yorbit ( [Segransan et al.|2()TT] ) by 
scrambling the RV data many times and counting the number 
of randomly generated peaks in the corresponding periodograms 
with an amplitude higher than the peak of the original set. Here, 
the 40-day signal of the original data was systematically higher 
than the peaks of 10000 scrambled sets. Longer period signals 
are also found around 120, 180, and 400 days. They are similar 
to the aliases of the one-year observational period of the spectral 
window. The middle panel of Fig. shows the periodogram of 
the RV constant master sampled at the same observing dates as 
HD 164595. We note that except for the 40-day triplet, other fea¬ 
tures are the same in the two periodograms. After correction of 
the systematic drift, that is, after subtracting the constant master 
RVs from the RVs of HD 164595, the only significant feature in 
the periodogram (FigJ^ bottom panel) is the 40-day triplet with 
a FAP lower than 10“ .This signal is then clearly not introduced 
by our correction, while other peaks were successfully reduced 
or removed. 

The multiplicity of the 40-day feature is due to the con¬ 
volution with the one-year observational period of the spectral 
window. The automatic de-aliasing algorithm CLEAN ( [Roberts | 
jet al.|198'7] ) implemented in Yorbit is another way to confirm the 
unicity of the physical signal in the feature by cleaning up the 
periodogram. While giving an important diagnostic, this method 
cannot be used to determinate which frequency is the physical 
one. This problem is particularly essential for the 40-day and 



2 . 5 . 10.0 20 . 50 . 100.0 200 . 500 . 

Period [days] 




2 . 5 . 10.0 20 . 50 . 100.0 200 . 500 . 


Period [days] 

Fig. 5. Upper panel : Generalized Lomb-Scargle (GLS) periodogram 
of HD 164595 RV before the constant master correction, with 10“^ 
FAP level (dashed line). Middle panel : GLS periodogram of the zero- 
point drift correction applied to HD 164595 dates of observation. Bottom 
panel : GLS periodogram of the corrected HD 164595 RV data, with 
10“^ FAP level (dashed line). 


45-day solutions, which are equally satisfying (with similar am¬ 
plitudes in the GLS periodogram and a similar ;^2 and residual 
dispersion for their respective fits). 

Dawson & Fabryd^ ( |201Q| ) proposed a method for solving 
this problem. The concept is to simulate a sinusoidal signal at the 
period and amplitude of each observational alias and to sample 
it at the date of observations. By comparing their periodograms, 
the one that significantly better matches the data periodogram 
corresponds to the true physical frequency. If the test is incon¬ 
clusive, it means that the noise prevents identification, and new 
observations are required. This was the case when we apply this 
method using only our 66 high-quality observations. We then 
used our complete data set of 75 measurements (see Fig.[^ and 
determined that the physical frequency is the 40-day peak. 

5.3. Distinguishing steiiar activity 

We studied the possibility that the 40-day signal might be due 
to star-induced activity and correspond to the rotational period. 
While we noted a slight increase in activity level with time, 
HD 164595 remained a relatively quiet star even in the last ob¬ 
serving season with a logR^j^ below -4.80. Figure [^ shows the 
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Fig. 6. GLS periodogram of the complete HD 164595 RV data (with all 
the 75 measurements) (solid black curve), the 40-day sinusoid (dashed 
blue curve) and the 45-day sinusoid (dotted and dashed red curve). The 
40-day sinusoid clearly is the best match to the RV data. 



Fig. 7. HD 164595 bisector Generalized Lomb-Scargle periodogram 
with 0.1 FAR level (dashed line). Long-period signals above this level 
are caused by a trend in the last season that is marginally correlated with 
the radial velocities. 


periodogram of the bisector span, which can be used to reflect 
stellar activity. No signal is present at 40 days or at its first har¬ 
monics. This also remains true with the FWHM and the logR'^^. 
There is no correlation between RV and these activity indica¬ 
tors either, except for the already mentioned marginal correla¬ 
tion between RV and the bisector, but only in the last season. 
Moreover, the rotational period of HD 164595 estimated from 


the log usin g the calibration published by Mamajek & Hil 
lenbrandl 120081) is 20 + 5 days, which is quite distinct trom the 


observed 40-day signal. We note that no 20-day signal appears in 
the RV periodogram, which means that the 40-day signal cannot 
be an observational alias. The HIPPARCOS photometry does not 
show any significant variations at a level of 13 mmag (with 147 
points), and with no hint of periodic signal in the periodogram. 
We also checked that the 40-day signal maintains the same phase 
and amplitude over the last two seasons when fitted indepen- 
dantly (the first season does not have enough data points to con¬ 
clude). We are therefore confident that the 40-day signal is not 
induced by activity and is planetary in origin. 


5.4. Radial-velocity analysis and orbital solution 

For the rest of the analysis and the fit of a Keplerian signal we 
only kept the high-quality data set (66 RV measurements) after 
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Table 4. Orbital and physical parameters of HD 164595b 


Parameter 

HD 164595b 

P [days] 

To [BID] 

e 

oj [deg] 
K[ms-i] 
msini [M©] 
a [AU] 

drift [m/s/yr] 

40.00 + 0.24 
56280 ±12 

0.088-^o'2 

145-16^ 

3.05 ± 0.41 
16.14 + 2.72 
0.23 

-2.34 + 0.44 

^meas 

66 

Data span [days] 

809 

o-(O-C) [ms-i] 

2.3 

Jitter [ms“^] 

1.8 


removing the nine low-quality measurements (low S/N, moon 
light contamination, anormal simultaneous thorium-argon flux, 
see Sect. 2). The RMS of the RV time series shown in Fig. 
(upper panel) is 3.73 ms“^ A slight long-term trend is visible. 
We used FASTIS to fit a Keplerian and a linear drift to our data 
without a priori. We ran 20 MCMC chains with 300 000 steps 
each, and discarded the solutions that did not converge to the 
40-day period (some converged to the observational aliases, in 
particular the 45-day period). 

Table reports the orbital parameters and their uncertain¬ 
ties for a Keplerian model with a linear trend, and the best fit is 
plotted in Fig. The fitted model of the planet has a 40- day 
period and a semi-amplitude of 3.05 + 0.41 ms“^ with an in- 
significa nt eccentricity. Assuming a mass of 0.99 ± 0.03 M© for 
the star ( |Porto de Meho et ar||2014| ), the minimum mass of the 
planet is 16.1 ± 2.7 M©. With a jitter of 1.8 ms“^ the reduced 
is equal to 1. The residuals after the fit present a dispersion 
of 2.3 ms“^ , and no significant feature appears in the corre¬ 
sponding periodogram. An interesting trend can be observed in 
the last season. It might correspond to either the activity-related 
drift that is correlated with the bisector span, or to a shift in the 
global trend and a subsequent incomplete correction. Trying to 
correct this trend on the data does not significantly change our 
reported solution. 

Another analysis using the genetic algorithm built into Yorbit 
( [Segransan et al.|20TT] ) yielded nearly identical parameter values 
and errors (P = 39.99 ± 0.19 days, K=2.96 ± 0.40 ms“^ and 
e=0.078±0.14). We note that the analysis of the uncorrected data 
does not change our results, but only increases the error bars of 
our orbital solution, the scatter of the residuals (from 2.2 to 2.9 
ms“^), and the corresponding jitter (from 1.8 to 2.7 ms“^). 


6. Discussion and conciusion 

We have reported the discovery of HD 164595b with the SOPHIE 
spectrograph. This source is a warm Neptune-like exoplanet with 
a minimum mass of 16.1 ± 2.7 M© orbiting a bright solar-analog 
star with a period of 40.00 ± 0.24 days with no significant eccen¬ 
tricity. An additional long-term RV drift is seen on our data. At 
this stage, we cannot conclude on its origin. It might be due to 
a long-period additional companion in the system, but it might 
likewise be the signature of a stellar magnetic cycle since we ob¬ 
serve an slight increase of the activity index over o ur observing 
span of 2.2 years. According to |Lovis et al.| ( [20ira| , a change of 
0.1 dex in the logR'^^ may introduce a change of up to 10 m s“^ 
in radial velocity. At this stage, we cannot completely exclude 
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Fig. 8. Upper panel : Radial velocities of HD 164595 with the best-fit or¬ 
bital solution (solid line). Error bars include the photon and instrumental 
noise and the stellar jitter. Bottom panel : Phase-folded radial velocities 
of HD164595 with the best-fit orbital solution (solid line). Error bars 
include the photon and instrumental noise and the stellar jitter. 


an under-correction of the zero-point drift either. Additional SO¬ 
PHIE measurements will help to conclude on the long-term drift. 

To validate and illustrate the capacity of SOPHIE to reach 
the 2 ms“^ precision thanks to an appropriate correction of sys¬ 
tematic instrumental effect, we refined the orbital parameters of 
the known planetary system HD 190360 based on SOPHIE and 
HIRES radial velocity measurements. With this precision on all 
timescales, SOPHIE can detect exoplanets with an RV semi¬ 
amplitude as low as 3 ms“^ Nevertheless, pushing the RV pre¬ 
cision down to 1 ms“^ will be crucial to efficiently cover the 
parameter space of low-mass exoplanets. Additional improve¬ 
ments are ongoing on SOPHIE, such as the thermal control of 
the spectrograph tripod and the implementation of a Fabry-Perot 
etalon within the calibration unit to improve the determination 
of the simultaneous drift. 

About 70 exoplanets with msin i smaller than 40 M® and an 
uncertainty on msin i smaller than 30% are known. They are dis¬ 
played in the mass - orbital period diagram in Fig. In that 
sample, no transiting low-mass planets wit h period greater than 
ten days are known except for Kepler-20c ( [Gautier et al.|2012) , 
which has a period of 10.9 days. Two-thirds of the 18 exoplan¬ 
ets with orbital periods longer than ten days are part of a multi¬ 


ple planetary system. Only three transiting planets with masses 
in the range 10-20 M® are known (Kepler-20c, GJ3470b, and 
HAT-P-26b). Their radii range from 3.1 to 6.3 R® , illustrating 
the diversity in density of this type of objects. This also under¬ 
lines the current lack of constraints on these objects. The transi¬ 
tion between super-Earths and Neptune-like planets and its de- 
pendance on the orbital period is poorly understood. The actual 
number of planets with precise density measurements is indeed 
far too low to cover the parameter space and constrain interior, 
formation, or evolution models for the low-mass planet popula¬ 
tion. The future transit search missions around bright stars TESS 
([Ricker et al. [2014[ ) and CHEOPS ( [Broeg et~ar[[2013t [Fortier 


[et al.|2Q14[ ) will address this issue from two complementary per¬ 
spectives. While TESS will substantially increase the number 
of transiting exoplanets suitable for high-precision RV measure¬ 
ment, CHEOPS will mesure precise radii of known exoplanets 
that have been discovered with radial velocities to derive precise 
densities. However, the low overall transit probability of the non¬ 
transiting planets of this sample (-300%) implies that CHEOPS 
will need a significant number of targets at various periods, re¬ 
quiring continuous efforts on low-mass RV surveys. 

In this context, HD 164595b is an interesting and potential 
target for the follow-up transit mission CHEOPS, although its 
transit probability is only 2%. HD 164595 matches all CHEOPS 
requirements in terms of coordinates, magnitude, and spectral 
type. The uncertainty of the transit window extrapolated to 2018 
is about six days at this stage, but long-term monitoring of this 
target with SOPHIE is planned to reduce this uncertainty. Fol¬ 
lowing the mass-radius relation given by [Marcy et al.[ ( [MT^ and 
assuming that sin /=1, the expected size of HD 164595b is around 
4.7 R®. Hence the expected transit probably has a depth close to 
1.9 mmag with a duration of 6.1 hours. We note, however, that 
the three transiting planets with similar masses as HD 164595b 
that were previously mentioned have a radius between 3.1 and 
6.3 R®. In case of transit, CHEOPS will be fully appropriate 
to derive the radius with a precision better than 10%. Finally, 
CHEOPS could also be used to detect transits of additional short- 
period low-mass companions that are unseen in radial velocities 
in the HD164595 system. 
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Table 5. Radial-velocity measurements, error bars and bissector for HD 164595. The dates with an asterix correspond to low quality data discarded 
for the analysis. 


BID 

RV[kms-M 

Uncertainty [ km s ^ ] 

Biss [kms 

56051.596 

2.1088 

0.0013 

-0.0248 

56058.576 

2.1094 

0.0013 

-0.0240 

56059.487 

2.1056 

0.0014 

-0.0282 

56061.524 

2.1112 

0.0018 

-0.0260 

56062.557 

2.1057 

0.0015 

-0.0285 

56072.524 

2.1118 

0.0013 

-0.0262 

56074.524 

2.1118 

0.0014 

-0.0283 

56076.503 

2.1117 

0.0013 

-0.0295 

56432.514 

2.1091 

0.0017 

-0.0270 

56434.581* 

2.1123 

0.0014 

-0.0270 

56436.544 

2.1077 

0.0018 

-0.0240 

56438.527 

2.1145 

0.0018 

-0.0238 

56443.518 

2.1100 

0.0016 

-0.0288 

56447.595 

2.1026 

0.0014 

-0.0253 

56449.414 

2.1022 

0.0014 

-0.0265 

56451.437 

2.1048 

0.0017 

-0.0273 

56454.408 

2.1032 

0.0015 

-0.0262 

56456.381 

2.1039 

0.0013 

-0.0290 

56462.520 

2.1048 

0.0016 

-0.0262 

56465.499 

2.1049 

0.0015 

-0.0277 

56468.582 

2.1067 

0.0017 

-0.0280 

56469.377 

2.1085 

0.0019 

-0.0290 

56470.421 

2.1093 

0.0015 

-0.0247 

56471.402 

2.1095 

0.0018 

-0.0238 

56486.475 

2.1099 

0.0015 

-0.0290 

56487.540 

2.1099 

0.0016 

-0.0278 

56489.510 

2.1075 

0.0014 

-0.0258 

56495.437 

2.1032 

0.0013 

-0.0280 

56497.549 

2.1034 

0.0014 

-0.0303 

56499.486 

2.1042 

0.0014 

-0.0267 

56501.441 

2.1059 

0.0013 

-0.0262 

56502.480 

2.1041 

0.0020 

-0.0327 

56517.400 

2.1114 

0.0014 

-0.0252 

56520.402* 

2.1097 

0.0013 

-0.0283 

56521.390 

2.1045 

0.0014 

-0.0270 

56524.431 

2.1102 

0.0019 

-0.0288 

56730.676* 

2.1143 

0.0013 

-0.0347 

56733.663* 

2.1124 

0.0017 

-0.0248 

56734.647 

2.1037 

0.0014 

-0.0280 

56753.586 

2.1137 

0.0014 

-0.0258 

56755.611 

2.1098 

0.0013 

-0.0222 

56756.578* 

2.1169 

0.0023 

-0.0287 

56758.605 

2.1112 

0.0013 

-0.0273 

56761.624 

2.1082 

0.0013 

-0.0240 

56777.522 

2.1037 

0.0015 

-0.0283 

56778.570 

2.1034 

0.0016 

-0.0200 

56782.561 

2.1047 

0.0012 

-0.0298 

56783.568 

2.1039 

0.0012 

-0.0313 

56786.585 

2.1081 

0.0014 

-0.0287 

56788.581 

2.1076 

0.0014 

-0.0267 

56792.607 

2.1072 

0.0014 

-0.0288 

56793.555 

2.1065 

0.0013 

-0.0295 

56794.603 

2.1078 

0.0012 

-0.0313 

56795.529 

2.1070 

0.0013 

-0.0342 

56815.447 

2.1043 

0.0011 

-0.0312 

56815.483 

2.1046 

0.0011 

-0.0323 

56816.560 

2.1007 

0.0012 

-0.0323 

56816.592 

2.1031 

0.0012 

-0.0327 

56817.421 

2.1020 

0.0011 

-0.0323 

56817.593 

2.1062 

0.0011 

-0.0282 
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56819.464 

2.1047 

0.0013 

-0.0322 

56819.564 

2.0997 

0.0012 

-0.0332 

56820.440 

2.0995 

0.0011 

-0.0310 

56820.514 

2.1013 

0.0011 

-0.0335 

56834.466 

2.1039 

0.0013 

-0.0295 

56838.411 

2.1085 

0.0015 

-0.0278 

56841.415 

2.1080 

0.0013 

-0.0322 

56844.519 

2.1062 

0.0012 

-0.0313 

56847.487 

2.1010 

0.0014 

-0.0320 

56850.546* 

2.0970 

0.0016 

-0.0363 

56854.401 

2.1044 

0.0012 

-0.0315 

56860.493 

2.0982 

0.0021 

-0.0422 

56900.431* 

2.0949 

0.0012 

-0.0310 

56901.312* 

2.1002 

0.0013 

-0.0317 

56902.312* 

2.0951 

0.0015 

-0.0360 
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Table 6. Radial-velocity measurements and error bars for HD 190360 (after correction). 


BID 

RV [kms-M 

Uncertainty [ km s ^ ] 

55756.500 

-45.2113 

0.0013 

55757.487 

-45.2097 

0.0012 

55758.507 

-45.2077 

0.0011 

55759.532 

-45.2057 

0.0012 

55760.522 

-45.2055 

0.0011 

55761.349 

-45.2096 

0.0012 

55762.493 

-45.2104 

0.0011 

55763.380 

-45.2151 

0.0012 

55765.622 

-45.2167 

0.0011 

55767.627 

-45.2197 

0.0014 

55769.479 

-45.2178 

0.0011 

55772.359 

-45.2144 

0.0011 

55785.487 

-45.2181 

0.0011 

55793.383 

-45.2056 

0.0011 

55796.467 

-45.2088 

0.0011 

55798.493 

-45.2107 

0.0012 

55813.428 

-45.2038 

0.0011 

55818.330 

-45.2123 

0.0010 

55835.417 

-45.2154 

0.0010 

55848.299 

-45.2048 

0.0010 

55908.236 

-45.2075 

0.0014 

56175.409 

-45.1906 

0.0012 

56200.333 

-45.1865 

0.0012 

56220.319 

-45.1856 

0.0011 

56517.451 

-45.1907 

0.0011 

56519.428 

-45.1901 

0.0011 

56521.443 

-45.1921 

0.0012 

56554.416 

-45.1969 

0.0011 

56557.378 

-45.1977 

0.0011 

56583.355 

-45.1921 

0.0012 

56613.268 

-45.1914 

0.0011 

56624.292 

-45.2050 

0.0011 

56625.297 

-45.2042 

0.0010 

56626.280 

-45.2035 

0.0011 

56628.290 

-45.1997 

0.0011 

56629.263 

-45.2015 

0.0011 

56630.238 

-45.1979 

0.0010 

56631.246 

-45.1935 

0.0010 

56638.266 

-45.2044 

0.0011 

56643.254 

-45.2092 

0.0010 

56644.242 

-45.2058 

0.0011 
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